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We study the statistical properties of population dynamics evolving in a realistic two-dimensional 
compressible turbulent velocity field. We show that the interplay between turbulent dynamics 
and population growth and saturation leads to quasi-localization and a remarkable reduction in the 
carrying capacity. The statistical properties of the population density are investigated and quantified 
via multifractal scaling analysis. We also investigate numerically the singular limit of negligibly small 
growth rates and delocalization of population ridges triggered by uniform advection. 
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For high nutrient concentration on a hard agar plate, 
the Fisher equation P can be a good description of 
the spreading of microorganisms such as bacteria at low 
Reynolds number Q. However, many microorganisms, 
such as those in the ocean, must find ways to thrive and 
prosper in high Reynolds number fluid environments. In 
presence of a turbulent advecting velocity field u(x, t). 
the Fisher equation reads 



dC_ 
'dt 



V • (uC) DV'^C + /iC(l - C), 



(1) 



where C (x, t) is a continuous variable describing the con- 
centration of micro-organisms, D is the diffusion coeffi- 
cient and ^ is the growth rate. As an example of "life 
at high Reynolds number" , we could take Eq. ([1]) to rep- 
resent the density of the marine cyanobacterium Syne- 
chococcus Q under conditions of abundant nutrients, so 
that /i ~ constant. 

As discussed in an advecting compressible turbu- 
lent flow leads to highly non-trivial dynamics. Although 
the results of Q were obtained only in one dimension 
using a synthetic advecting flow from a shell model of 
turbulence, two striking effects were observed: the con- 
centration field C(x, t) is strongly localized near tran- 
sient but long-lived sinks of the turbulent flows for small 
enough growth rate /x; in the same limit, the space-time 
average concentration (denoted in the following as carry- 
ing capacity) becomes much smaller than its maximum 
value 1. Both effects are relevant in biological applica- 
tions [HI. 

In this Letter, we present new numerical results for 
more realistic two dimensional turbulent flows. We as- 
sume that the microorganism concentration field C(x, i), 
whose dynamics is described by Eq. ([1]) , lies on a planar 
surface of constant height in a three dimensional fully 
developed turbulent flow with periodic boundary condi- 
tions. Such a system could be a rough approximation 



to photosynthetic microorganisms that actively control 
their bouyancy to manitain a fixed depth below the sur- 
face of a turbulent fluid [1]. As a consequence of this 
choice, the flow fleld in the two dimensional slice becomes 
compressible Q • We consider here a turbulent advecting 
field u(x, t) described by the Navier-Stokes equations, 
and nondimcnsionalizc time by the Kolmogorov dissi- 
pation time-scale = (z^/e)^/^ and space by the Kol- 
mogorov length-scale rj = [v^ /e)^/^^ where e is the mean 
rate of energy dissipation and v is the kinematic viscosity. 
The non-dimensional numbers charecterizing the evolu- 
tion of the scalar field C(x, t) are then the Schmidt num- 
ber Sc = v/D and the non-dimensional time ^r^. A 
particularly interesting regime arises when the doubling 
time Tg = is somewhere in the middle of the in- 
ertial range of eddy turnover times (r^ = r/6rU, where 
SrU is the typical velocity difference across length scale 
r) that characterize the turbulence. Although the under- 
lying turbulent energy cascade is somewhat different , 
this situation arises for oceanic cyanobacteria and phyto- 
plankton, who double in 8 — 12 hours, in a medium with 
eddy turnover times varying from minutes to months . 

The main results of our investigation are the follow- 
ing: we confirm the qualitative behaviour found in [3| 
for a two-dimensional population, under more realistic 
turbulent flow conditions. We also investigate the limit 
/i — > cxD and discuss the singular limit of /i — > validating 
the physical picture proposed in Our understanding 
of the limit /xr^ ^ 1 may be helpful in future investi- 
gations, as explicit computations for /i > can be very 
demanding. In addition, we quantify the statistical prop- 
erties of the concentration field and investigate the effect 
of an uniform convective background flow field. 

We conducted a three dimensional direct numerical 
simulation (DNS) of homogeneous, isotropic turbulence 
at two different resolutions (128'^ and 512'^ collocation 
points) in a cubic box of length L = 27r. The Taylor mi- 




2 




FIG. 1: (Color online) (Left panel) Pseudocolor plot of concentration field. The bright yellow regions indicate regions of 
high concentration (C > 0.1) and the black regions indicate regions of low concentration. (Right panel) Pseudocolor plot of 
[C(x, to)/(0.1 + C(x, fo))(tanh(— V • u))]. The dark red regions indicate negative divergence and large concentration whereas 
dark blue regions indicate positive divergence and large concentration. Plot are made at identical time to (after the steady 
state has been reached) on a slice z = const obtained from our 512^ numerical simulations of Eq. ^ for /ir^ = 0.0045 and 
Schmidt number Sc = 5.12. Note that microorganisms cluster near regions of compression (V • u < 0), as is evident from the 
high density of red regions. 



croscalc Reynolds number for the full 3D simulation 
was Rex = 75 and 180, respectively, the dimcnsionless 
viscosities were v = 0.01 and i/ = 2.05 • 10""^, the total 
energy dissipation rate was around e ~ 1 in both cases. 
For the analysis of the Fisher equation we focused only 
on the time evolution of a particular 2D slab taken out 
of the full three dimensional velocity field and evolved a 
concentration field C (x, t) constrained to lie on this plane 
only. This is a particularly efficient way of producing a 
compressible 2d velocity field in order to mimic the flow 
at the surface of oceans. Note that our velocity field has 
nothing to do with 2d turbulence and has the structures, 
correlations and spectra of a bidimensional cut of a fully 
3d turbulent flow. A typical plot of the 2d concentration 
field and the concentration field conditioned by the corre- 
sponding velocity divergence field (taken at time t ~ 86, 
Re\ = 180) in this plane is shown in Fig. [1] {Sc = 5.12). 
The Fisher equation was stepped forward using a second- 
order Adams-Bashforth scheme. The spatial derivatives 
in the diffusion operator are discretized using a central, 
second-order, finite-difference method. As the underlying 
flow field is compressible, sharp gradients in the concen- 
tration field can form during time-evolution. In order 
to capture these sharp fronts we use a Kurganov-Tadmor 
scheme for the advection of the scalar field by the velocity 
field 0. 

The concentration C (x, t) is highly peaked in small 



areas, resembling one dimensional filaments (see Figure 
[1] and supplemental movie). When the microorganisms 
grow faster than the turnover times of a significant frac- 
tion of the turbulent eddies, C(x, t) grows in a quasi- 
static compressible velocity field and accumulates near 
the regions of compression, leading to filaments [l^ • The 
geometry of the concentration field suggests that C(x, t) 
is different from zero on a set of fractal dimension dp 
much smaller than 2. A box counting analysis of the 
fractal dimension of C(x, t) supports this view and pro- 
vides evidence that dp = 1.0±0.15. 

A biologically important quantity is the spatially aver- 
aged carrying capacity or the density of biological mass 
in the system, 

Z{t) = ^ J dxdyCi^;t), (2) 

and in particular its time average in the statistical steady 
state as a function of the growth rate fi, (Z)^. With- 
out turbulence (Z)^ = 1 for any fi. When turbulence 
is acting, in the limit /i — > cxd we expect the carry- 
ing capacity attains its maximum value (Z)^^oo = 1, 
because when the characteristic time Tg becomes much 
smaller than the Kolmogorov dissipation time t^, the ve- 
locity field is a relatively small perturbation on the rapid 
growth of the microorganisms. Indeed, consider a per- 
turbation expansion of C(x, t) in terms of Tg. We define 
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C(x, t) = ^ T^C„(x, i), and substitute in Eq. ([T]), 

n— 0,...,oo 

where the functions C„ (x, t) are the coefficients of the 
expansion. Upon assuming a steady state and cohecting 
the terms up to 0{Tg) we find, after some algebra, 

{Z), « 1 - (r2/i^)(/(V • u)Mx) + 0(r3). 

The hmit /i ^ can be investigated by noting that for 
ji — Q, Eq. ((ij reduces to the Fokker-Planck equation 
describing the probabihty distribution V{x,t) to find a 
Lagrangian particle subject to a time varying force field 
u(x,t): 

dV 

— + V • (u-p) = ZJV^P. (3) 
at 

Upon defining F = ((V • u)^)^/^ as the r.m.s value 
of the velocity divergence, following Q we expect a 
crossover in the behavior of {Z)f^ for fi < T. In the limit 

— ?• 0, we expect: 

To understand Eq. note first that for small ^ the 
statistical properties of C should be close to those of P. 
Thus we can assume that, in a statistical sense, C(x, t) w 
L'^{Z)fj,V{x,t). Averaging Eq. ([T]) in space and time 
leads to (C)^ - (C^)^ = 0, which is equivalent to Eq. (g]). 
Eq. (|4]) is crucial, because it allows us to predict for 
small /i from the knowledge of the well-studied statistical 
properties of Lagrangian tracer particles without growth 
in compressible turbulent fiows. We have therefore tested 
both Eq. (j4]) and the limit — > oo against our numerical 
simulations. In Fig. [2] we show the behavior of (Z)^ for 
the numerical simulations discussed in this Letter. The 
horizontal line represents the value 1/((7'^)L^) obtained 
by solving Eq. ([3]) for the same velocity field and /i = 0. 
The insert shows a similar result for a one dimensional 
compressible flow Q. For our numerical simulations we 
observe, for /ir,, > Ft^ « 0.23 the carrying capacity (Z)^ 
becomes close to its maximum value 1. The limit /.( — > 
requires some care: the effect of turbulence is relevant for 
Tg longer than the Kolmogorov dissipation time scale r^. 
We take the limit /i — !• at fixed system size L. When 
Tg ^ tl ~ (L^/e)^/'^, the large scale correlation time, 
the population is effectively frozen on all turbulent time 
scales, and Eq. ([4]) should apply. 

The limit /i — > can be investigated more precisely 
as follows: according to known results on Lagrangian 
particles in compressible turbulent flows, V should have 
a multifractal structure in the inviscid limit 1/^-0 0, 
If our assumption leading to Eq. (|4]) is correct, 
C (x, t) must show multifractal behavior in the same limit 
with multifractal exponents similar to those of V. 

We perform a multifractal analysis of the concentration 
field C{x,y,t) with /i > by considering the average 
quantity Cfj,{r,t) = Jg,. C{x,y,t)dxdy where B{r) is 



a square box of size r. Then the quantities {Cfj,{r)P) are 
expected to be scaling functions of r, i.e. {C^{r)P) ^ 
j-O'ip) ^ where a{p) is a non linear function of p. 

In Figure [2] we show the quantity a{p) for /z = 0,0.1, 
and 1 for < p < 4 extracted from power law fits to 
{C{ry) over ~ 1.5 decades. In the inset we show (C'(r)^) 
for fi = 0.01. Although our dynamic range is limited, the 
scaling description seems to work with smoothly vary- 
ing exponents a{p). Even more important, the statistical 
properties of (7^ (r) seems to converge to the case /Lt = as 

— > 0. In the same figure we also show, for comparison, 
a similar analysis performed for the energy dissipation 
field (black line); see for a detailed description. 

Our multifractal analysis suggests a relation between 
the quasi localization length ^ and the carrying capacity 
{Z)f^,. The quasi localization length ^ can be considered 
as the smallest scale below which one should not observe 
fiuctuations of C(x, t). In the limit /_i — > 0, we can define 
the quasi localization length ^ as: 

e = ^""'^ (5) 

~ ((V7')2)- ^ ^ 

We expect ^ to be of the same order of the width of 
the narrow filaments in Fig. [1] To compute (Z)^ as a 
function of ^, we observe that it is reasonable to assume 
(752 (x,t)) ~ {C^{r = ~ Using Eq. © we 

obtain (Z)^ - ^-"(2). Qn the left side of Figure [3] we 
show {Z)f^ as a function of ^ (obtained by using ([5])) for 
/i = 0.01 by varying the diffusivitiy D. Reducing the 
diffusivity D shrinks the localization length ^ and (2')^ 
becomes smaller. Dimensional analysis applied to Eq. 
[3] and suggests that ^ oc ijfiv j t)^!^ . More generally 
we expect that ^(-D) is monotonically increasing with D. 
From Figure 13] (left side), a reasonable power law behav- 
ior is observed with a scaling exponent 0.46 ± 0.03 very 
close to the predicted behavior — a(2) = 0.47 obtained 
from Fig. [2] (right side inset). 

Note that Fig. [3] (left side) reveals a strong depen- 
dence of (Z)^ on ^ and hence on the microbial diffu- 
sion constant. D in turn depends on the ability of ma- 
rine microorganisms to swim. Approximately 1/3 of the 
open ocean isolates of S'ynechococcus can propel them- 
selves along their micron-sized long axis at velocities of 
'-^ 2bpLml sec [isj . Upon assuming a random direction 
change every '--^ 20 — 30 body lengths, the effective diffu- 
sion constant entering Eq. ([T]) can be enhanced 1000— fold 
relative to D for passive organisms. The extensive energy 
investment required for swimming in a turbulent ocean 
becomes more understandable in light of the increased 
carrying capacity associated with say, a ~ 30— fold in- 
crease in ^. Some marine microorganisms may have 
evolved swimming in order to mitigate the overcrowding 
associated with compressible turbulent advection. 

Finally, we discuss bacterial populations subject to 
both turbulence and uniform drift because of, e.g., the 
ability to swim in a particular direction. In this case, we 
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FIG. 2: (Color online) Left panel: Behavior of the carrying capacity {Z)fj,T^ as a function of fiTr, from 128^ (red dots) and 
512^ (blue squares) numerical simulations with Sc = 1. Note that for /ir^ < 0.001, the carrying capacity approaches the 
limit 1/((P^)L'') ^ 0.16 ± 0.02 (blue line) predicted by Eq. (g}. In the inset we show similar results for one dimensional 
compressible turbulent flows in [J]. Right panel: The anomalous exponents a{p) computed by the multifractal analysis of 
the concentration field C(x,t) for different values of and Sc = 1 for our 512^ numerical simulation. Note that for — >■ 
the multifractal exponents approach the statistical properties of the field V described by Eq. (|3J. The black line shows the 
multifractal properties of the energy dissipation rate e. In the inset we show the scaling behavior of (C'(r-)^) for /i = 0.1. 



can decompose the velocity field into zero mean turbu- 



lence fluctuations plus a constant drift velocity uq [14 | 
along e.g. the x-direction. In presence of a mean drift 
velocity Eq. ([T|) becomes: 



uoe^)C] = DV'^C + mC(1 - C) (6) 
where Bx is the unit vector along the cc-direction. Note 
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FIG. 3: (Color online) Left: Log-log plot of {Z) as a function 
of the localization length ^ defined in Eq. ([SJ for ua — 0. The 
blue line is the fit to the power-law. The slope is consistent 
with the prediction (Z) ^ ^""(2) (jisgygsed in the text. The 
numerical simulations are done for ^ = 0.01 and different 
values of D from D = 0.05 to D = 0.001. Right: Plot of 
{Z) as function of a super-imposed uniform velocity uo for 
fj, = 0.01 (red bullets) and /i = 0.1 (green triangles) with 
D = 0.015. 



that a mean drift (to follow nutrient gradients, say) 
breaks the Galilean invariance as the concentration C 
is advected by ug, while turbulent fluctuations u remain 
fixed. In Fig.|3]we show the variation of carrying capacity 
versus uq for two different values of /i and fixed diffusiv- 
ity D = 0.015. We find that for uq ^ Wrms (the root- 
mean-square turbulent velocity) the carrying capacity Z 
saturates to a value equal to the value of Z in absence 
of Uq i.e., quasilocalization by compressible turbulence 
dominates the dynamics. For uq > Urms the drift veloc- 
ity delocalizes the bacterial density eventually causing 
Z ^ 1, as was also found in d = 1. 

Wc have shown that a realistic model for two dimen- 
sional compressible turbulence predicts reduced microor- 
ganism carrying capacities, similar to those found in a 
highly simplified Id model Simulations at two ele- 
vated Reynolds number show that results are robust and 
in agreement when properly normalized. The limit of 
large growth rates was addressed analytically, and the 
statistics maps smoothly onto known results for con- 
served densities advected by compressible turbulence. Fi- 
nally we studied the effect of a preferred swimming di- 
rection on the carrying capacity. 
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